Hidden Order in the Cuprates 
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We propose that the enigmatic pseudogap phase of cuprate superconductors is characterized by a 
hidden broken symmetry of d x i_ y i-type. The transition to this state is rounded by disorder, but in 
the limit that the disorder is made sufficiently small, the pseudogap crossover should reveal itself to 
be such a transition. The ordered state breaks time-reversal, translational, and rotational symme- 
tries, but it is invariant under the combination of any two. We discuss these ideas in the context 
of ten specific experimental properties of the cuprates, and make several predictions, including the 
existence of an as-yet undetected metal-metal transition under the superconducting dome. 
PACS numbers: 71.10.Hf, 71.27.+a, 74.72.-h, 71.10.Pm 



I. INTRODUCTION 

In this paper we argue that much of the strange phe- 
nomenology of the cuprate superconductors may be sim- 
ply explained as the disorder-frustrated development of 
a new order parameter. There are a number of potential 
candidates for this order, but the one we favor on phe- 
nomenological grounds is orbital antiferromagnetism 
or d-density wave (DDW) order [||^] , which is character- 
ized by a local order parameter which distills the univer- 
sal physics underlying the staggered flux state ||-^| di- 
vorced from the uncontrolled approximations associated 
with the gauge theory formalism. The essence of our 
idea is that the pseudogap || observed in underdoped 
cuprates is an actual gap in the one-particle excitation 
spectrum at the wavevector (7r,0) and symmetry-related 
points of the Brillouin zone associated with the devel- 
opment of this new order. It is "pseudo" in experiment 
only because of extreme sensitivity to sample imperfec- 
tion caused by proximity to the phase transition. More- 
over, the DDW couples weakly to common experimental 
probes, and is thus difficult to detect. 

Our proposal has much in common with theoretical 
ideas already in the literature [[s|-|l9|], and borrows heav- 
ily from them. For example, Wen and Lee have pro- 
posed staggered currents that fluctuate but do not or- 
der Q . Varma has proposed currents which alternate in 
the unit cell but do not break translational symmetry 
@. Emery and Kivelson 0@ and Caprara et al. Q 
have proposed states with broken symmetries of different 
kinds. Our strategy for constructing a theory and con- 
fronting experiments differs from most others in deem- 
phasizing modeling of the "strange metal" behavior and 
focusing on order, low-temperature phenomenology, and 
material imperfection - all issues with sharp experimen- 
tal dichotomies amenable to falsification. DDW order 
can be detected if it is present. If it is not present, the 
proposal is disproved. 
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1. Top: Phase diagram constructed from experiments 
on YBa2Cu307_«. Tjv is the Neel transition, T c is the su- 
perconducting transition, T* is the pseudogap crossover [27], 
To is the location of the maximum in the uniform suscepti- 
bilty, and Tchargc is a charge-ordering line recently reported 
by Mook et al. [39]. Middle: Values of \z\ (solid), y (dots) 
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(solid) minimizing the free energy of Eq. 



(1) with A = 1, 7 = —0.8, and the linear functions a and a' 
shown on the bottom. The parameter p is hole doping. 
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FIG. 2. Left: Arrangement of bond currents in the DDW 
state. Right: Brillouin zone of the Cu-O plane. The dots 
show the half-filling Fermi surface, as well as the Brillouin 
zone boundary after the DDW state has formed. The circles 
are the Fermi surface of the DDW state at finite doping. 



II. COMPETING ORDER 

Order-parameter competition has always been a nat- 
ural candidate for explaining why the superconducting 
transition temperature T c first grows and then retreats 
as doping is reduced. Let us consider the generic zero- 
temperature Ginzburg-Landau free energy 

F = X(y 2 + \z\ 2 ) 2 +1 y 2 \z\ 2 -ay 2 -a'\z\ 2 , (1) 

describing the development of order parameters y and z 
in the case that low-order mixing is forbidden by sym- 
metry. In Fig. lb we plot the values of y and z that 
minimize F for the case of A = 1 and 7 = —0.8 as a func- 
tion of the abstract tuning parameter p. The variables a 
and a' are the simple linear functions of p shown in Fig. 
lc. One sees that z develops at p = 0.3, y develops at 
p = 0.2, and that 0.1 < p < 0.2 is a coexistence region 
in which the growth of y suppresses and eventually elim- 
inates z. Thus if we imagine z to be the magnitude of 
the order parameter for d-wave superconductivity and p 
to be doping, then we can understand the onset, growth, 
saturation, and eventual destruction of superconductiv- 
ity with reduced doping as an effect of a monotonically 
strengthening d-wave pairing interaction, as opposed to 
one that first strengthens and then weakens. The un- 
derdoped side of the superconducting dome is then fun- 
damentally different from the overdoped side in that the 
superfluid density is suppressed there by the development 
of a second order parameter y. 



III. D-DENSITY WAVE 

Let us now consider the order parameter 



(2) 



where /(k) = cos(k x ) — cos(k y ). If /(k) were replaced by 
a function with s-wave symmetry, y would simply be the 
order parameter of a charge-density wave (CDW) - hence, 
we call this state a d x 2_ y 2 density wave state (DDW) ||. 
For the particular case of Q — (tt,tt), which we think 



most relevant to the cuprates, the equivalence of Q and 
— Q enforced by the underlying band structure requires 
the sum to be imaginary. Thus this state necessarily 
breaks parity and time-reversal symmetry (i.e. exhibits 
magnetism) , as well as translation by one lattice spacing 
and rotation by 7r/2. It is, however, symmetric under the 
combination of any two of these operations. The order 
parameter is equivalent to the array of bond currents 
illustrated in Fig. 2. 

The excitation spectrum of the DDW at very low en- 
ergies is generic and consists of conventional fermionic 
particles and holes in a band structure like that of the 
d-wave superconductor with which it competes. Intro- 
ducing a mean-field ansatz fl (cf. Eq. ^) we obtain the 
1-body Hamiltonian 



(3) 



where ek = — 2t[cos(k x )+cos(k y )} and Ak = yV[cos(k x ) — 
cos(fcy)], and V is a coupling constant in the micro- 
scopic Hamiltonian. Microscopic Hamiltonians with 
short-range repulsion and superexchange are favorable 
for such order j^,^] but are even more favorable for an 
antifcrromagnetic state. However, correlated hopping 
terms tend to tip the balance in favor of DDW order 
||. Since the ordering occurs at Q = (ir,7r), it is most 
favorable at half-filling or low doping. The corresponding 
band structure is 



E k = ± 



| A DDW| 



(4) 



At half-filling there are gapless quasiparticles only at the 
nodal points k = (±7r/2, ±7r/2). At finite doping, Fermi 
pockets are opened, as shown in Fig. 2. While the DDW 
state is semimetallic at half-filling, it is a conventional 
metal (with 2-d localization prevented by interlayer tun- 
neling) with a disconnected Fermi surface at dopings 
other than half-filling. It is possible for the DDW to 
discommensurate, thereby opening a full gap, as occurs 
with a traditional spin density wave, but this is not auto- 
matic because the remaining Fermi surface is not nested. 
Some related density- wave states are discussed in the Ap- 
pendix. 

The excitation spectrum at high energies is not generic. 
There is no reason for the quasiparticlc at (tt, 0) to have 
integrity, particularly if the system is near the continuous 
quantum phase transition at p — 0.2 in Fig. 1. This is 
a Fermi-surface reconnection, at which the Hall conduc- 
tance jumps, a van Hove singularity develops at (w, 0), 
and quasiparticles scatter violently even at low energy 
scales lOfillQ. 



IV. D-WAVE SUPERCONDUCTIVITY 

The Heisenberg exchange nominally responsible for 
DDW order also tends to favor d-wave superconductivity. 
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This is the underlying reason the band structures of the 
two are so similar, and why the competition of these two 
kinds of order is natural. If we allow the superconduct- 
ing bond expectation value < CjfCjtj, >= ±z to develop, 
where the sign is positive on x bonds and negative on y 
bonds, the Hartree-Fock Hamiltonian becomes 

H' HF = Hhf +J ^2 ±{z 4 ; ct T + z* CjfCki) , (5) 
<jk> 

and the corresponding superconducting quasiparticle dis- 
persion relation becomes 



H 



[(4 



| A DDW|2)l/2 ±/i ]2 + |AD SC | S 



(6) 



where A[? sc = zJ[cos(k x ) — cos(fc y )] and /i is the chemi- 
cal potential. Thus not only does this kind of interaction 
stabilize both kinds of order, it allows the two order pa- 
rameters to evolve continuously into each other without 
collapsing the quasiparticle gap at the zone face. This 
allows us to use the ground state expectation value of 
TL and similar Hamiltonians as a sensible model for the 
energy functional F, i.e. one that does not throw away 
important low-energy excitations. 

This calculation illustrates an important feature of the 
mixed state that the superfluid density is not fixed by 
sum rules on the underlying Fermi surface but is rather 
determined by the balance between the DDW and DSC 
order parameters. This is because the superfluid is pri- 
marily a condensate of Cooper pairs drawn from the 
gapped region near (it, 0) rather than the residual Fermi 
surface near (tt/2,it/2). This effect is not difficult to un- 
derstand if the DDW order parameter is small, for then 
the semimetallic state with Fermi points - or, away from 
half-filling, the conventional metallic state with a small, 
disconnected Fermi surface - is not significantly differ- 
ent from the parent metal with a full Fermi surface at 
the energy scales relevant to superconductivity. As the 
DDW order parameter becomes large, however, we have 
more and more the case of a powerful attractive force 
exciting electrons and holes virtually into the insulat- 
ing part of the band structure and then binding these 
into superfluid. The result is a condensate fraction that 
falls precipitously as the DDW order parameter grows. 
An insulating ground state (or, in this case, nearly in- 
sulating, since there is a small disconnected Fermi sur- 
face) that becomes a superfluid without first becoming 
a metal is unusual in solids, but perhaps not in nature, 
for this is the central idea behind Higgs condensation in 
electroweak theory. 



V. S-WAVE COMPETITION 

The competition between DDW and DSC has a sim- 
ple analogue in the s-wave cas e p3| that is particularly 
instructive because it is exact 1241. The Hubbard model 



<jk> a j 



(7) 



has the special property at half-filling that replacing the 
fermion operator on a lattice site j, Cjj_, by (— l) J 'cL (a 
unitary transformation at half-filling) reverses the sign of 
U. When U > this model has a ground state which is 
an ordered antiferromagnet characterized by the expec- 
tation values 



<Sf> 

<5 i > 
<S]> 



< c jt c n 



> 



1 < C tf C 3l - C 

< c h c n- c n c n> 



(8) 



When U < the ground state is thus a degenerate 
mixture of s-wave superconductivity and checkerboard 
charge order characterized by the expectation values 



<Re(A J )> 
<Im(A J )> 



* <C 3l C 3l 



i-iy <c] T cj t + c] l c J i> 



(9) 



Both kinds of order occur simultaneously, are equivalent 
energetically, and may be rotated into each other by anal- 
ogy with spin rotation of an antiferromagnet. More pre- 
cisely, this system lies at a quantum phase transition be- 
tween the two kinds of order and can be made to acquire 
one, the other, or a mixture of the two by means of an 
arbitrarily small perturbation, exactly the way the pa- 
rameters 7 and a — a' in Eq. (1) break the rotational 
invariance of F. 

The Hartree-Fock solution, which is only approximate, 
also has this symmetry. Allowing the expectation values 
y = (— I y < c^Cji — ct^Cji > /2 and z =< cj^Cji > we 
obtain for the Hartree-Fock Hamiltonian 

Hhf = -tJ2 E c Lc fe . + (-l) J y(c] T c jT 

<jk> jo j 



(10) 



and for the corresponding quasiparticle dispersion rela- 
tion 



E k = ± 



J[(el + |ACDW|2)l/2 ± ^2 + | A SSC|2 ( ( n ) 



where A CDW = yU and A ssc = zU, are the charge den- 
sity wave and s-wave superconducting gaps, respectively. 
This is exactly the same as Eq. @ with s-wave quantities 
substituted for d-wave ones. Thus, as in the d-wave case, 
the superconducting order parameter may, at half-filling 
(fi = 0) be rotated continuously from pure superconduc- 
tivity to pure checkerboard charge order without closing 
the quasiparticle gap. In this case, however, the rotation 
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also leaves the ground state energy invariant, and is an 
exact symmetry [ fj"6|f24| . 

This calculation illustrates the important feature of 
charge order that it competes easily and naturally with 
s-wave superconductivity but not with <i-wave. This is 
because it is an s-wave condensate, per Eq. (||). 

VI. PSEUDOGAP 

A large number of experimental properties of the 
cuprates are consistent with the presence of DDW order 
in underdoped samples. 

A. Gap Evolution 

The ci-wave superconducting gap in the electron spec- 
tral function evolves continuously with underdoping into 
the cZ-wave-like pseudogap without collapsing. In the top 
of Fig. 3 we reproduce point-contact tunneling measure- 
ments on underdoped YBCO of Renner et al. |2f| show- 
ing the excessive size of the tunneling gap and its persis- 
tence above the superconducting T c , both of which are 
characteristic of underdoped cuprates. Identification of 
this feature with the d-wave gap follows from its evolu- 
tion out of the simpler BCS-like gap found in overdoped 
materials and its rough compatibility with the magnitude 
of T c . However, its persistence above T c is not consistent 
with a traditional BCS gap, for this should disappear 
at T c , as occurs in overdoped samples, on quite general 
grounds. That this gap has the correct angular depen- 
dence is shown in the middle of Fig. 3, where we repro- 
duce angle-resolved photoemission spectra from under- 
doped Bi2Sr2CaCu20g+5 at two different points in the 
Brillouin zone reported by Norman et al. pq| . The up- 
per trace, taken from near the zone face at (71", 0), shows 
a large gap that persists well above T c , whereas the lower 
trace, taken from near the node at (7r/2,7r/2), shows a 
smaller gap that is destroyed at T c . This angular depen- 
dence is also seen in the bottom of Fig. 3, where we repro- 
duce the retreat of the photoemission "leading edge" as a 
function of position on the weak-coupling Fermi surface 
reported by Harris et al. p7{ . The d- wave-like character 
of the gap is clear, as is its persistence at the zone face 
above T c for even slightly underdoped samples. Thus it 
appears that the pseudogap and the superconducting gap 
have identical functional forms and evolve continuously 
into each other as the doping is reduced, just as expected 
from order-parameter rotation. 

Energetic competition as the cause of this rotation is 
suggested by the similarity between the superconducting 
and pseudogap energy scales. It may be seen in the bot- 
tom of Fig. 3 that the maximum "leading-edge" gapis 
30 meV while the retreat caused by heating above T c is 
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FIG. 3. Top: Point-contact tunneling spectrum of under- 
doped Bi2212 by Renner et al. [24] at temperatures (top to 
bottom) 4.2 K, 63 K, 81 K, 89 K, 109 K, and 151 K. The 
dotted line shows the spectrum at T c — 85 K. Middle: An- 
gle-resolved photoemission from Norman et al. [25] on under- 
doped Bi2212 with T c — 75 K at the two points in the Brillouin 
zone shown in the insets. The data were symmetrized to re- 
move the Fermi function. The temperatures are from top to 
bottom: 65 K, 85 K, 110 K. Bottom: Leading-edge gap from 
Harris et al. [26]. Left panel: T c = 78 K Dy-BSCCO, at tem- 
peratures 13 K (diamonds), 100 K (crosses), 150 K (squares). 
Right panel: T c = 46 K Dy-BSCCO; the symbols are the 
same as in the left panel. 

between 5 meV and 10 meV, depending on doping. The 
pseudogap scale k^T* ~ 30 meV is also identified in 
a number of other measurements S, notably neutron 
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FIG. 4. Top: Specific heat 7 in mj/mole K 2 for various 
doping levels of Bi2.i5Sri.85CaCu20s+<5 of Tallon and Loram 
[31]. The y-intercepts of the curves increase with doping so 
that the lower curves correspond to underdoped crystals while 
the upper curves correspond to overdoped crystals; the fourth 
highest curve corresponds to optimal doping. Bottom: Spe- 
cific heat jump [f{T c ) - 7(120 K) ] of above samples ver- 
sus doping. The doping level is determined from T c and the 
semiempirical relation relation between this and doping p in 
holes per Cu site shown in plusses. 

scattering ^8|, NMR |^{|, electronic Raman scattering 
pOf , and optical reflectivity |pl[ . 



B. Superfluid Density 

Rapid collapse of superfluid density below optimal 
doping is seen in many experiments fl32| . The zero- 
temperature penetration depth, for example, grows 
rapidly in the pseudogap regime and correlates with 
the suppression of T c with underdoping, yet saturates 
at overdoping in a way reminiscent of a traditional 
BCS superconductor Q. In Fig. 4 we reproduce the 
heat capacity measurements on Bi 2 isSri.ssCaC^Os+s 
recently reported by Tallon and Loram p2] ]. Above a 
hole concentration of about p — 0.19 per Cu the spe- 
cific heat jump at the superconducting transition varies 
weakly with p, as one would expect if the material were 
an ordinary metal undergoing a transition to BCS super- 
conductivity. At p — 0.19, however, there is an abrupt 
transition and a rapid decrease of this height with under- 
doping, as though all or part of the Fermi surface were 
being destroyed by the removal of holes. As a result of 
this, there are fewer low-energy excitations remaining to 
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FIG. 5. Inelastic neutron scattering at k = (n,n,n) re- 
ported by Bourges [33] for YBCO at various levels of under- 
doping. The closed and open circles correspond to T = 5K 
and T = 100K, respectively. Higher-resolution experiments 
[30] have now shown that the shoulder on the low-energy side 
below T c is actually split incommensurate peaks. 

be affected by the superconducting transition. Hence, the 
specific heat jump at the transition, A7, is reduced. All 
of this behavior is compatible with Fig. 1 if the phase 
transition is at p = 0.2, where the order parameter y 
begins to develop, is associated with the onset of DDW 
order and the consequent continuous opening of a gap at 
(tt, 0) in the quasiparticle spectrum. 



C. Spin Susceptibility 

There is evidence that spin ordering - and thus pre- 
sumably stripe ordering - has not taken place at optimal 
doping in YBCO, but only occurs at much lower dop- 
ing levels. In Fig. 5 we reproduce the inelastic neutron 
measurements [p4| for optimally-doped and underdoped 
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E 

FIG. 6. Im Xq( tJ ) as defined by Eq. 13 for the case of 
H = 0, yJ/(2t + xJ) = 0.1 and U/U c = 0.0,0.3,0.6,0.9, where 
U c is the critical value of U. 



and broadens, just as occurs with decreased doping in 
Fig. 5. This width is due to efficient decay of the exci- 
ton into nodal quasiparticle pairs. At a slightly higher 
value of U the continuum evolves into a divergence at 
uj = associated with onset of spin order. Note that the 
DDW and DSC order parameters in this calculation are 
effectively interchangeable. Since A£? DW = — Aj^^ and 
A£ sc = — for q = (7r,7r), the coherence factor is 

unity and unchanged close to the Fermi energy whether 
or not both gaps, or only one of them, are present. For 
an s-wave gap the corresponding coherence factor would 
have been zero. 



YBCO at a momentum transfer of (tt,tt,tt). These ex- 
periments show that the 41 meV resonance, which dis- 
appears above T c and is presumably associated with the 
superconductivity, evolves continuously with underdop- 
ing into the magnetic fluctuation spectrum of the or- 
dered antiferromagnet. Thus, we interpret the piling 
up of low-frequency spectral weight in the experiment at 
low doping as signaling the approach of magnetic order, 
and conversely of showing that magnetic order is neither 
present nor imminent at the onset of DDW order. The 
spin-fluctuation spectrum in the superconducting region 
remains fully gapped and has no low-energy structure 
of any kind. The resonance continues to be destroyed 
by elevated temperature, but the requisite temperature 
grows with underdoping even as T c is evolving to zero. 
In this way an excitation manifestly associated with the 
superconductivity at optimal doping transforms into an 
excitation irrelevant to superconductivity |3jJ. 

This effect is simply understood as a triplet exciton 
p6| that vanishes at elevated temperature because the 
quasiparticle gap required for it to be well-defined van- 
ishes. This is quantified in Fig. 6, where we plot the 
imaginary part of 



1 + U XqM 



(12) 



where 



k+q 



(w + ii]) 2 - (E k + £ k+q ) 



\DDW \ DDW , /\DSC ADSC\ 

^k+q + ^k ^k+q I Qg\ 



at q = (tt, 7r) for various values of U. This is a crude 
ladder sum in which x q ( w ) represents the susceptibility 
of the ideal BCS superconductor characterized by -Ek, 
A£ DW , and A£ sc per Eq. (§), while U represents a 
coulomb interaction added to push this system toward 
spin antiferromagnetism. One sees that as U is increased 
the sharp resonance in the spectrum decreases in energy 



D. High-Field Transport 

Stripes and antiferromagnetic order are naturally asso- 
ciated with the insulating behavior of the cuprates seen 
near half-filling j37| . In a conventional doped band insu- 
lator, insulation is caused by impurities, which trap carri- 
ers and prevent them from moving. The system becomes 
a metal when it is doped sufficiently that the impurity 
orbitals touch. One of the most significant characteristics 
of the cuprates is that they continue to insulate to phe- 
nomenally high dopings, typically 5% or 1 hole for every 
20 Cu atoms. It is very difficult to understand how an in- 
sulator with an energy gap less than that of the common 
semiconductor GaAs should still insulate at these high 
dopings through impurity trapping solely. But develop- 
ment of antiferromagnetic order with antiphase domain 
walls, which then trap carriers and pin, is easy to un- 
derstand, physically sensible, and supported experimen- 
tally by the simultaneous occurrence in these materials 
of discommensurated magnetic Bragg peaks and X-ray 
satellites at exactly half their momentum displacements 
[B8[ . Thus our view is that charge ordering (which would 
have an order parameter of the form (^J), but with an 
/(k) which has s-wave symmetry and, in all likelihood, 
incommensurate Q) impedes conduction, rather than fa- 
cilitating it juj, and moreover is characteristic of the 
insulating state. 

The issue of coexistence of superconductivity with 
stripes and antiferromagnetism, and potential causative 
relations among them, is still highly controversial and 
a matter of experimental study ]3l| . There is, how- 
ever, increasing evidence that the coexistence found 
in I^-xSr^CuO^Nd |Q is anomalous and that the 
cuprates with the highest values of T c have charge or- 
dering only at the low-doping edge of the supercon- 
ducting dome. The recent neutron scattering from 
YBa 2 Cu0 3 7 _ x reported by Mook et al. [Q find the 
charge-ordering line shown in Fig. 1 and no static anti- 
ferromagnetism anywhere in the superconducting region. 
This is consistent with the the high-field transport exper- 
iment on Bi2Sr2- a; La ;r Cu06+<5 recently reported by Ono 
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FIG. 7. High-field transport experiment of Ono et al. [40] 
on Bi2Sr2- a; La a: Cu06+i • Top: Resistivity in mf2-cm versus 
temperature for samples with La concentrations (top to bot- 
tom) x= 0.84, 0.76, 0.73, and 0.39 with (circles) and without 
(solid lines) an applied magnetic field of 60 T. Bottom: The 
temperature of the minimum value of p (diamonds) plotted 
versus hole concentration. The superconducting T c (dots) is 
plotted for reference. 

et al. |4lf reproduced in Fig. 7, which finds a metal- 
insulator transition at essentially the same doping as the 
charge ordering line of E3 ] when the superconductivity 
is crushed by a large magnetic field. The phenomenol- 
ogy of this transition is qualitatively similar to that ob- 
served previously in La2_ K Sr x Cu04 fi^] except that it 
occurs near the edge of the dome rather than near opti- 
mal doping. This is important, for Castellani, Di Castro, 
and Grilli ji^] were led by this observation to propose 
that the strange-metal behavior of the cuprates might be 
quantum criticality associated with the charge-ordering 
transition. These more recent experiments suggest that 
it is instead quantum criticality associated with the de- 
velopment of DDW order. LSCO is unique among the 
high-T c cuprates in having a low transition temperature, 
a strong tendency to stripe-order near 1/8 doping, and 
an extreme sensitivity to Nd doping pS[ ], all of which 
suggest mechanical weakness of the crystal structure. 

The large-field experiment also reveals another impor- 
tant aspect of the cuprates, namely the lack of evidence 
for strange-metal behavior in the zero-temperature nor- 
mal state. It may be seen in the top of Fig. 7 that 
the resistivities on the metallic side of the transition be- 
come constant at low temperatures and that they evolve 
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FIG. 8. c-axis optical conductivity of YBa2Cu30e.7 in 
fl~ 1 cm~ 1 reported by Homes et al. [43] for temperatures 
(top to bottom) T = 250 K, 150 K, 110 K, 70K, and 10 K . 
Large phonon contributions have been subtracted out. 

continuously across the transition into the linear- T re- 
sistivity characteristic of the high-temperature normal 
state of the cuprates. The resistivity at the transition 
is also about 200 /if2-cm, a typical saturation resistivity 
in strong-scattering metals. Both of these properties are 
consistent with the zero-temperature normal state being 
a conventional metal. They do not prove this, but they 
make the argument for a non-Fermi-liquid phase more 
difficult, as linear-T resistivity is one of its key signa- 
tures. 

Thus on the basis of these experiments we predict that 
in large magnetic fields there should be a second zero- 
temperature phase transition near p = 0.19 associated 
with the onset of DDW order. At this transition the 
system should remain a conventional metal but violently 
change the topology of its Fermi surface. This transition 
should be plainly visible in all transport measurements 
and should be characterized by powerful critical scatter- 




E. c-Axis Conductivity 

DDW formation provides a simple explanation for 
the perplexing semiconducting c-axis resistivity in many 
cuprates. In Fig. 8 we reproduce the optical conductiv- 
ity measurements on YBa2Cu306.7 of Homes et al. |44|] 
showing the steady reduction of the oscillator strength 
below 40 meV beginning at a temperature far above the 
superconducting T c . That this reduction does not con- 
serve the f-sum rule locally - which any mean-field the- 
ory, including that of the DDW, does - is interesting but 
not necessarily significant, as a mean-field description is 
obligated to be quantitative only at arbitrarily small en- 
ergy scales. Band structure studies |l5| of these materi- 
als have shown that the c-axis tunneling matrix element 
is largest at (71", 0) and symmetry-related points - pre- 
cisely at the points where the DDW gap is large. (The 
functional form is roughly t± ~ (cosfc^ — cosk y ) 2 Ea| .) 
Thus the opening of the DDW gap suppresses the c-axis 
transport because the remaining Fermi surface does not 
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conduct efficiently in the c-direction due to small tun- 
neling matrix element. The above matrix element holds 
for simple tetragonal materials (Hgl201, T11201, etc.). 
For body centered tetragonal materials (LSCO, T12201, 
Bi2212, etc.), the maxima of t± are shifted towards the 
zone center, and the effect of opening the DDW gap at 
(tt, 0) and symmetry-related points is weaker. 



VII. ORBITAL MAGNETISM 

The distinguishing characteristic of DDW order is the 
magnetic field it makes. Since the possibility of sponta- 
neous breaking of time-reversal and parity in the cuprates 
was first proposed in the late 1980s there have been a 
number of attempts to detect such fields, most of which 
have reported null results . However there has always 
been confusion about the size of the effects one would 
expect, and there have always been mysterious magnetic 
signals in the cuprates, including a recent report of spin 
antiferromagnetism coexisting with superconductivity in 
a sample of superoxygenated La2Cu04+j, with y=0.12 
and T c = 42 K This fundamentally conflicts with 

a previous report of no magnetism in Lai.ssSro.isCuCU 
p6| . We feel that the magentic experiments are so con- 
tradictory that they can at present neither rule out nor 
confirm the presence of DDW order. 

We estimate the magnetic field at the center of a pla- 
quette associated with DDW order to be between 1 and 
30 gauss jlj. The bond currents of Fig. 2 are roughly 
eA DDW /S, where A DDW is the maximum DDW gap. If 
we take this to be 30 meV, we find bond currents of 
about 7/iA The large uncertainty in the corresponding 
field strength is due mainly to uncertainty in the cur- 
rent path. One can reasonably consider models ranging 
from Cu sites connected by 1 A "wires" to split current 
carried between adjacent O atoms. 

Let us now briefly review the current experimental sit- 
uation relevant to direct detection of DDW magnetism. 



A. Neutron Scattering 

DDW order is, in principle, visible in magnetic neu- 
tron scattering. Unfortunately the signals are quite small 
compared with those from ordered spins and easily over- 
whelmed by them. The ratio of the staggered magnetic 
field associated with DDW fields to that nominally pro- 
duced by an ordered array of spins is 



B 



DDW 



-Bafm 



3 A DDW \ fmer" 



her 



eh 



^A DDW , (14) 



or about 0.06, with r = 4 A taken for the bond length. 
Effective magnetic moments of this size are just barely 
detectable in the cuprates fl39}f47| ] . 



It is also unfortunate that the doping levels at which 
DDW order should be well developed lie close to the spin- 
glass regime p§] where the system crosses over between 
Neel and superconducting order. The spin glass is char- 
acterized by slightly incommensurated short-range anti- 
ferromagnetism with strongly suppressed scattering in- 
tensities along one orthorhombic axis - behavior consis- 
tent with unpaired Cu spins pointing in the plane |47| ] 
and inconsistent with DDW magnetism. However, nu- 
merous incursions of this magnetism into the supercon- 
ducting phase have been reported, in one case deeply |39| , 
and this has always been difficult to understand from the 
point of view of traditional magnetic models. It implic- 
itly raises the question of whether there might be two 
kinds of antiferromagnetism in the cuprates - one, due to 
spins, which is incompatible with superconductivity and 
one, due to DDW, which is fully compatible with it and 
associated with pseudogap formation. Spin-orbit cou- 
pling would then mix these and conceivably make them 
evolve into each other with increased doping. 



B. X-ray Scattering 

DDW order cannot be seen in X-ray scattering. The 
DDW order parameter is odd under time-reversal while 
atomic displacements are even, so there is no first-order 
coupling between them, and Bragg scattering through 
circular birefringence from the valence electrons is too 
weak. For a 10 KeV X-ray of frequency u> the Bragg 
intensity is down by the factor (fi B B DDW jhwf ~ 10~ 16 
from the Bragg intensity of valence electrons - already 
small compared with the signal from the core electrons. 
The absence of an X-ray signal is a key characteristic 
DDW order distinguishing it experimentally from stripes. 



C. Magnetic Resonance 

The static magnetic field of ideal DDW order cannot 
be seen directly through NMR of Cu or O nuclei in ideal 
CuO planes, as these lie at centers of symmetry where the 
DDW magnetic field is zero. However, magnetic fluc- 
tuations associated with the onset of DDW order or a 
glassy state of a disorder-frustrated DDW could be seen 
by NMR, although it would be difficult to distinguish 
from antiferromagnetic spin fluctuations for the reasons 
stated above. Also, DDW order can in principle be seen 
in NMR of ions out of the Cu-0 planes, such as Y, Ba, 
La, or Sr. It has long been established that there are 
unusual magnetic signals below T c in all the cuprates, 
but attempts to quantify these have been plagued by 
the inherent model-sensitivity of NMR analysis. Tallon 
and Loram Q have recently argued using the ratio of 
63 Cu and 17 spin-lattice relaxation rates analyzed with 
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FIG. 9. Cuprate phase diagram inferred by Niedermayer et 
al. [47] from ^tSR measurements. The solid and open symbols 
refer to Yi_ I Ca I Ba2Cu306 and I^-xSr^CuCU, respectively. 
Tjv is the Neel temperature, T c is the superconducting tran- 
sition temperature, Tf is a secondary magnetic ordering that 
occurs on top of spin antiferromagnetism, and T g is the glass 
transition temperature. The dephasing time for the latter is 
approximately 0.1 /J.S. 

the model of Millis, Monien, and Pines |50| that short- 
range antiferromagnetic fluctuations develop in the pseu- 
dogap regime with a functional dependence on p tracking 
roughly the value of y in Fig. 1. This analysis is not per- 
suasive evidence for DDW order. 

D. Muon Spin Resonance 

Muon spin resonance has consistently found evidence 
for magnetism in the superconducting state of the 
cuprates for dopings less than p = 0.1. In Fig. 9 
we reproduce the phase diagram of Niedermeyer et al. 
|f3lf showing boundaries of distinct magnetic behav- 
iors observed in powders of both I^-^Sr^CuC^ and 
Yi_ 2; Ca a ;Ba2Cu3O6.02. They report a "spin freezing" 
transition (Tf) below the Neel transition, and a spin- 
glass transition (T g ) that cuts in to the superconduct- 
ing dome. Below this transition, and at doping levels 
as high as p — 0.09, the muons depolarize in about 0.1 
/xS. In the case of LSCO, the measurements extended 
into the range of the 1/8 anomaly at p = 0.12, beyond 
which the spin-glass neutron signal tends to disappear 
and where no magnetism was found in previous /iSR 
measurements [|4(| . However the fact that both cuprates 
behave similarly, and that the spin-glass temperature in 
Yi_ a; Ca a ;Ba2Cu306 is substantially higher, suggests that 
this behavior is characteristic of the cuprates as a class. 
Also, the way the spin-glass line ends has always been 
confusing. 

In a recent paper Panagopoulos et al. [|52| have 
reported anomalous long-time magnetic fluctuations 
at temperatures just above the glass transition in 
La2_ a; Sr a ;Cu04 powders. The scale of these is compa- 
rable to T c and has a functional dependence with dop- 
ing identical to that of the parameter y in Fig. 1 - i.e. 
decreasing with doping and vanishing at p — 0.17. Thus 
they argue that the spin-glass line actually ends here, not 



at 1/8. The observation of the same effect in a different 
cuprate, which seems likely in light of Fig. 9, would sug- 
gest an intrinsic magnetic signal developing at the onset 
of DDW order. 

VIII. DISORDER AND CROSSOVER 

The muon phenomenology suggests an answer to a 
question plaguing the idea of competing order in the 
cuprates, namely why there is no evidence for a genuine 
phase transition at the pseudogap temperature T*, the 
alleged phase boundary for onset of DDW, and also why 
previous searches for magnetism at optimal doping have 
found sample-dependent or null results. It is simply that 
the DDW order is corrupted by disorder and transformed 
into the spin-glass transition at the lower temperature T g 
p3| . In very dirty samples it is lowered so much as to be 
effectively destroyed. 

There has been controversy over how much intrinsic 
disorder cuprates possess since they were discovered. The 
essence of the problem is that the most sensitive tests of 
disorder - transport and the degradation of superconduc- 
tivity - are corrupted by the non-Fermi-liquid behavior 
of the normal state |2(|] evidenced by resistivities which 
exceed the Ioffe-Regel limit of 100 /xf2-cm at T c (c/. Fig. 
7) and increase with temperature from there. However, 
using criteria less dependent on the theory of metals, the 
case for chronic disorder is easier to make: All cuprates 
lose oxygen easily in arbitrary amounts. All of them have 
spin-glass phases at low dopings [Q. All of them have 
magnetic scattering in the superconducting regime that 
is sample-dependent, difficult to reproduce, and difficult 
to quantify [^(J. All of them have anomalous widths in 
Cu and O NMR and NQR Q. Thus our view is that 
all cuprates made thus far have been significantly disor- 
dered, even ones showing evidence to the contrary such as 
narrow superconducting transition widths. This view is 
supported by new scanning tunneling microscope exper- 
iments on atomically perfect cleaves of optimally-doped 
BSCCO that find inhomogeneities in the tunneling den- 
sity of states on the scale of 20 A |35) . 

We note that the disorder need not occur within the 
CuC>2 planes to have a strong effect on the electronic 
properties. The non-superconducting cuprates which 
were studied in the late 1980s differ from the supercon- 
ducting ones only in the elements that sit between planes. 
Substituting Hg between the planes raises T c substan- 
tially. Substituting Nd causes stripes (38). 

The DDW transition is in the same universality class 
as the random-bond Ising model. The DDW order pa- 
rameter breaks translational and rotational symmetries, 
and thus couples to disorder as in a model with a random 
uniaxial anisotropy. From the Imry-Ma argument |53[ as 
adapted to the random anisotropy case [ p6[ , we know 
that this symmetry-breaking transition will be spoiled 
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by the random distribution of impurities. Thus, in the 
presence of disorder, time-reversal is the only true sym- 
metry that can be spontaneously broken by the DDW 
state. The universality class is then that of a Zi symme- 
try preserved by the impurities. 

The phase diagram of the random-bond Ising model 
depends critically on the disorder strength. Weak dis- 
order is an irrelevant perturbation and can be ignored 
at the finite-temperature transition to a state with bro- 
ken time reversal symmetry. Such a state has a non- 
vanishing expectation value for the staggered orbital 
magnetization (thereby breaking the disorder- averaged 
translational symmetry). If the disorder is strong, on 
the other hand, and the interlayer coupling is finite, then 
there can be a finite-temperature transition in the same 
universality class as the three-dimensional (3D) Ising 
spin-glass transition. Due to the weakness of the mag- 
netic coupling between the planes, the spin-glass tran- 
sition temperature estimated from the two-dimensional 
spin-glass susceptibility f^l , \ 2 S g ~ T~ 7 , 7 ~ 5.3, is 
small. Such a transition would not be possible if wc 
could neglect the coupling between the planes, as the 
lower critical dimension of the Ising spin glass is known 
to be greater than two. 

In contrast to this, the finite temperature transition 
to DSC remains sharp in the presence of disorder - al- 
though T c may be degraded. This is because disorder 
does not couple linearly to the order parameter as a ran- 
dom field, and because a superconducting transition in 
two dimensions is possible. The ultimate 3-dimensional 
transition driven by the coupling between the layers is 
robust. This can be further understood by invoking the 
Harris criterion (5^] assuming that the transition is in the 
3D-XY universality class. The criterion states that weak 
disorder is an irrelevant perturbation to the pure system 
if the specific heat exponent is negative, which is indeed 
the case for the 3D-XY model. When the disorder is so 
strong that kpl ~ 1, where Uf is the Fermi wave vector 
and I is the mean free path, a superconductor-insulator 
transition will take place, and the Harris criterion will no 
longer apply. 

The potential presence in this system of a disorder- 
sensitive, purely electronic, phase transition involving a 
Fermi-surface reconnection raises the disturbing possibil- 
ity that many experiments in this field may be measur- 
ing corrupted critical properties of the DDW transition 
rather than the properties of new states of matter. The 
notorious non-Fermi-liquid behavior of the normal state, 
for example, appears to evolve at the lowest temperatures 
and in a strong magnetic field into behavior of a tradi- 
tional metal. A possible explanation of this is that the 
high-temperature behavior is characteristic of a quantum 
critical |5q] region associated with a nearby critical point. 



IX. SUMMARY 

In summary we find that most of the strange behavior 
of the cuprate superconductors is consistent evidence for 
the simultaneous occurrence of d-wave superconductiv- 
ity and bond antiferromagnetism. On the basis of this 
we predict that the spin-glass transition temperature ob- 
served in muon spin resonance will climb to the pseu- 
dogap temperature T* as the sample quality improves, 
that the onset of this effect with doping coincides per- 
fectly with the loss of superfluid density at p = 0.19, and 
that this transition will be found to be a metal-metal 
transition involving a Fermi surface reconnection, not a 
transition to stripe order pl|- p^j37| , |38[ | , when magnetic 
fields sufficiently intense are available to crush the super- 
conductivity at optimal doping. 
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APPENDIX: RELATED DENSITY WAVES 

There are two related unconventional density wave or- 
der parameters potentially relevant to the cuprates ||. 
The first is that the frustration of the singlet DDW order 
parameter can lead to incommensurate ordering, in anal- 
ogy with the Ferrell-Fulde-Larkin-Ovchinnikov |5^] state 
in superconductivity as nesting is destroyed. As in the 
superconducting case, this will take place for sufficiently 
strong frustration and at sufficiently low temperatures. 
Note that in this case the order parameter is allowed 
to couple with lattice displacements and can therefore be 
seen in X-ray scattering. When the order parameter is in- 
commensurate, it will no longer have pure d x 2_ y 2 symme- 
try, but will mix in p-wave terms. For Q = (ir/a, n/a) + q 
with |q| small, the order parameter will take the form of 
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equation (0), with /(k) = cos^) — cos(fc y ) replaced by [15] 



/(k) = (1 + -fea) [cos{k x a] 



cos(k y a)] 



- ^q x a sm(k x a) - ^q y a sin(k y a) 



(15) 



The second interesting order parameter is the triplet 
version of the DDW This is defined by 



y 



= J2 /( k ) S ^ < C k+Q, s c k, S ' > 



(16) 



where a is a Pauli spin matrix. If /(k) were chosen to 
be a function of s-wave symmetry, this would be a con- 
ventional spin-density wave. The order in this case is 
chracterized by broken time-reversal, translational, and 
rotational invariances. The combination of any two of 
time-reversal, a translation by one lattice spacing, or a 
rotation by n/2 is preserved, however. In addition, spin- 
rotational symmetry is also broken, which leads to gap- 
less spin-1 excitations. The triplet DDW corresponds 
to an alternating pattern of spin currents analogous to 
charge currents of Fig. 2. Presently, the phenomenology 
of high temperature does not seem to be consistent with 
the choice of the triplet DDW as the competing order 
parameter. 
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